Order from Disorder in the Two- Dimensional Kondo-Necklace 
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We analyze the effects of site-dilution disorder on the thermodynamic properties of the two- 
dimensional Kondo necklace using finite-temperature stochastic series expansion. Results will be 
discussed for the dependence on dilution concentration, temperature, and Kondo exchange-coupling 
strength of the uniform susceptibility, the staggered structure factor, and the Chakravarty-Halperin- 
Nelson ratio. Dilution is shown to induce effective free-spin clusters in the gapped phase of the clean 
system with a low-temperature Curie constant renormalized below 1/4. Furthermore, dilution is 
demonstrated to generate antiferromagnetic order in the quantum disordered phase of the clean 
system, i.e. order-from-disorder. In turn, the quantum critical point of the clean system, separating 
an antiferromagnetic from a paramagnetic dimerized state at a critical Kondo exchange-coupling 
strength J c w 1.4 is suppressed. Finally, speculations on a renormalized classical behavior in the 
dilution induced ordered phase are stated. 

PACS numbers: 75.10.Jm, 05.70.Jk, 75.40. Cx, 75.40. Ms 



Quantum critical points (QCPs), i.e. zero temperature 
phase transitions as a function of some control parameter 
are likely to be at the core of unconventional finite tem- 
perature behavior of many novel materials^. Quantum 
antiferromagnets (AFMs) with an intrinsic spin dimeriza- 
tion and weak inter-dimer exchange like K(Tl)CuCl3^ 
or BaCuSi2 0g£, are of particular interest here, since they 
allow for switching between quantum disordered spin- 
gapped phases and states with magnetic long-range order 
(LRO). QCPs have been induced in these materials both, 
by tuning the inter-dimer exchange via pressur o^iTiS , as 
well as by applying external magnetic field o 9 ' 10 ? 11 . 

Combining the physics of quantum critical spin sys- 
tems with that of disorder is an open issue. Site dilution 
with non-magnetic impurities has been observed to in- 
duce LRO in the spin gapped phases of several dimerized 
quantum AFMa 12 i 13 i 14 i 15 . In spin ladders and dimerized 
spin chains a picture of weakly interacting 'defect mo- 
ments' has emerged with enhanced AFM correlations in 
the vicinity of the non-magnetic sites^iiiSii^. In di- 
mensions D > 2 and at finite impurity concentration this 
may trigger LRO in spin-gapped systems on bipartite lat- 
tices and remove a QCR2&. In the context of the cuprate 
parent compounds^i, site dilution of the 2D Heisenberg 
AFM has been shown to lead to a QCP at percolation 
where LRO is suppresse d 22 ' 2 '? . In addition to site dilution 
several other forms of disorder are of interest. Early on 
random exchange has been investigated in a variety of ID 
and 2D Heisenberg and Ising AFMs where it can lead to 
random singlet fix-point a 24 ' 2,1 ?' 2 ^ and Griffiths phases^. 
Very recently 'dimer dilution' has been studied in bi- 
layer Heisenberg models yielding multicritical points at 
percolatio n 2 fii 2 Mft . 

In this work we focus on random site dilution in the 
SU{2) symmetric 2D spin-1/2 Kondo-necklace (SKN) 
which is shown in fig. 1 
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= j ^ Sp; • Sp„, + J ^ Spi ■ Sj 
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Figure 1: 2D Kondo-necklace. Spin-1/2 moments coupled by 
j and J and located at vertices in 'conduction-plane' P and on 
'Kondo-sites' in plane I at end of dash-dotted lines, 'o' labels 
randomly diluted 'Kondo-sites' with no magnetic moment. 



with j = 1 hereafter. This model can be viewed as the 
strong coupling limit at large Coulomb correlations in the 
conduction electron band of the Kondo-lattice model at 
half filling2i. In turn, the two layers P and I in fig. ^ will 
be loosely referred to as 'conduction-band' and 'Kondo- 
site' layer hereafter. The clean limit of the 2D SKN has 
been studied at finite temperatures recentl y'? 1 "'? 2 . The 
model exhibits a QCP between AFM-LRO and a dimer 
phase at J c ~ 1.4. It shows temperature-scaling and 
a Chakravarty-Halperin-Nelson ratio consistent with the 
O (3) non-linear sigma model. The aim of this work is to 
shed light onto the effect of randomly substituting non- 
magnetic ions into the Kondo-layer, as shown in fig. 
In that case the site index I in in the second sum of eqn. 

runs only over those sites which are occupied by a spin 
in layer I. 

Our analysis is based on the stochastic series expan- 
sion (SSE) with loop-updates introduced by Sandvik 
and Syljuase nffftift 4 . Averaging over disorder configura- 
tions is performed by rapid thermalization using the 
temperature-halving scheme proposed by Sandvih^. We 
refer the reader to the latter three references for details. 
In the following, all of the SSE results reported for a finite 
dilution concentration comprise of an average over 1000 
disorder configurations on a system of TV = 24 x 24 x 2 
sites. The latter is justified by the rather weak finite size 
dependence observed beyond this for the clean system^!. 
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Figure 2: (a), (b): Uniform susceptibility \u at J = 2 > J c 
vs. temperature for c=0.03, 0.1 and 0.2 and 1/1024 < T < 1 
and for c = and 0.05 < T < 1 on a lin-lin (a) and a log- 
log scale (b). (c): T\u product at c / as compared with 
'free spin' approximation Txu,fs from eqn. 10, for identical 
parameters as in (a), (d): T\u product at c 7^ as compared 
to Curie behavior renormalized by 1/3. Statistical errors in all 
panels either less than the solid-circle marker size or depicted 
by error bars. Legends label plots from top to bottom. 



To begin we discuss the uniform susceptibility 

Xu = (m 2 ) /T , (2) 

where T is the temperature and m = J2n=p+i 1 ^ni/^° ls 
the total spin ^-component with N° being the total num- 
ber of sites occupied by a spin. The defect concentration 
is c = 2Nf/N , where Nf is the number of empty site 
in layer /, i.e. N° = (2 — c) N/2. In fig. [2 we summa- 
rize several aspects of Xu versus temperature at J > J c 
which refers to the gapped phase of the clean system. 
The figure contrasts the clean system for 0.05 < T < 1 
against three impurity concentrations c = 0.03, 0.1, and 
0.2 for 1/1024 < T < 1. First, fig. Efa) clearly demon- 
strates that the non-magnetic defects induce low-energy 
magnetic density of states in the spin gap. Second, the 
log-log plot of fig. EJb) shows, that the susceptibility of 
these states is similar to that of 'free' spins, i.e. it follows 
almost a Curie law in the temperature range below the 
spin gap. 

To further elucidate the free spin behavior, we start 
from the strong coupling limit, i.e. J — > 00, of the Kondo- 
necklace. In that case and in the clean system the ground 
state is a product of singlets on the dash-dotted bonds of 
fig. H Removing a 'Kondo' spin leaves a free spin in the 
'conduction' plane. Denoting the uniform susceptibility 
per occupied site by Xu.fs, with the index fs referring to 
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Figure 3: Uniform susceptibility Xu vs - temperature at c = 
0.03 for J = 0.8 (dashed-double-dotted), 1.4 (dashed), and 
2.4 (solid), corresponding to the AFM-LRO, the critical and 
the gapped regime at c = 0. Statistical errors are depicted by 
vertical bars. 



the free spins, one has 



Xu,fs 



2(l~c) 
(2-c) 



Xu (c = 0) + 



(2 
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4T 



(3) 



The first term is the contribution from the singlets and 
the second term stems from the free spins. To assess 
the relevance of this free spin picture - other than in 
the strong coupling limit - we may use eqn. (pj with 
Xu (c = 0) taken from the SSE at c = but for J other 
than 00. In fig. Efc) Txujs from this approximation 
is compared with Txu from eqn.JSJ at J = 2. First, at 
low temperatures Txu is smaller than Txujs- Second, 
while Txuj's saturates on an energy scale set by the spin 
gap, Txu continues to decrease as a function of T and 
seem to saturate at far lower temperatures. This is a 
clear indication of low-energy correlations which develop 
between the impurity induced degrees of freedom. This 
also clarifies the need for SSE results at far lower values 
of T at finite disorder as compared to the clean case. 

For site diluted spin-ladders Sigrist and Furusaki have 
proposed^ that below a crossover temperature, and due 
to long-range, higher order exchange the free spins will 
form clusters of a size set by the inverse temperature. 
These clusters lead to a Curie-type of susceptibility per 
defect, however with a Curie constant reduced from 1/4 
to 1/ (3 • 4) . The argument of refill is not restricted to 
ladders. To consider this, fig. G{d) shows the Curie 
contribution normalized to the defect concentration on 
a log-log scale. For c > 0.1 saturation of Txu can 
be anticipated for T < 1/500. Consistently with refill 
liniT^o [Txu (2 — c) /c] < 1/4. Moreover, while a renor- 
malization by a factor of exactly 1/3 can not be read off 
from this figure, the reduction is very close to this value. 
The kink in observable at elevated temperatures and 
small concentrations in this log-log plot is related to the 
opening of the spin gap. 

Fig. 03 clarifies the difference between the gapped and 



3 



10 



.o 



; (a)c = 0.03 

Mi 

iflli 

- * i 
'iii 


: (b)c = 0.1 

■ — ■ Curie 
— • Curie/3 

- M 

Tjlli 


] (C)C = 0.2 

— - J = 1.6 

— 1.8 
- — 2.0 

— 2.2 

f~~ 2.4 

l r 


Hi 






J i i 







0.1 

0.001 0.01 0.1 



001 0.01 0.1 0.001 0.01 0.1 
T 



Figure 4: Txu vs. T at c = 0.03 (a), 0.1 (b), and 0.2 (c) for 
various values of J > J c . Dashed straight lines correspond to 
(renormalized) Curie constant. Legends label plots according 
to order listed between dashed and dashed-dotted. Statistical 
errors are depicted by vertical bars. 
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Figure 5: Low-temperature staggered magnetization 35 (Q) 
versus J for various c = 0, 0.03, 0.1, and 0.2 (a) and versus 
c for various J = 0.8 . . . 2.4 (b). Statistical errors are either 
less than the solid-circle marker size or are depicted by vertical 
bars. Legends label plots from top to bottom. 



the ordered phase as well as the impact of finite size gaps. 
It contrast to the gapped phase, site dilution of the LRO 
phase is not expected to induce quasi free moments be- 
cause of the AFM correlations between all the occupied 
sites^i. On finite systems however, a size-dependend spin 
gap will remain also in the AFM-LRO phase. Disorder 
will eventually induce states in this gap, leading to an 
increase of the uniform susceptibility, however at a tem- 
perature scale much lower than and unrelated to J. This 
is consistent with the results shown for small concentra- 
tion in fig. ED In the gapped phase, at J=2.4 an almost 
straight-line behavior can be observed on a log-log scale 
below an energy set by the intrinsic spin gap. In the LRO 
phase at J=0.8 such behavior can only be anticipated at 
temperatures at least two orders of magnitude lower. 

In fig. |3|we summarize T\u versus T for c = 0.03, 0.1, 
and 0.2 over a finite range of coupling constants J corre- 
sponding to the gapped regime of the clean system. First, 
while low-temperature convergence to a constant T\u is 
observed for c = 0.1 and 0.2, fig. QJallows for an approxi- 
mate extrapolation to T = at c = 0.03 only. Second, it 
is tempting to claim that at fixed c the low temperature 
Curie constant depends only weakly on J. Third, the low 
temperature Curie constant is less than 1/4 and close to 
1/ (3-4) for all cases depicted. However, the figure sug- 
gests that the Curie constant is a decreasing function of 
c and cannot be described by a single renormalization 
factor at all c. This is consistent with the fact, that at 
c = 1 the SKN is identical to the 2D AFM Heisenberg 
model, for which lim-r^o T\u = 0. 

Now we turn to the longitudinal staggered structure 
factor, i.e. the order parameter for AFM-LRO 



S (Q) 



(4) 



Kl ex P (»Q ■ r l N ° is the dag- 
gered magnetization with Q = (71", tt, 7t). All results dis- 



where tUq 



cussed correspond to values of T such that low temper- 
ature saturation of S (Q) has been reached. Fig. Efa) 
shows the squared staggered moment Mq = 35* (Q) ver- 
sus J for the clean case and for c=0.03, 0.1 and 0.2. In 
the clean case Mq is finite below the critical value of 
,/ = J c and, apart from finite size effects, drops to zero 
for J > J c . I.e. J c corresponds to the QCR21 with AFM- 
LRO for J < J c and a spin gapped state for J > J c . 
Generating SSE data of smaller value of Mq for J ~ J c 
at c = requires 'fine-tuning' of J which we refrain from. 

The main result of fig. EI a) is contained in the change 
of Mq versus J upon doping. First, for all concentrations 
investigated the QCP disappears. Second, a finite stag- 
gered moment can be observed at all values of J and in 
particular also in the formerly gapped phase of the clean 
system. This implies that doping by non magnetic impu- 
rities induces AFM-LRO in the quantum disorder phase 
of the SKN, i.e. order-from- disorder. It is tempting to 
speculate, that this behavior is true for all c. 

We may also discuss these results from a different point 
of view, i.e. by considering Mq versus the impurity con- 
centration for various exchange coupling constants, as in 
fig. GJb). This shows Mq to increase at small c, both 
below and above J c . At J c the induced AFM order is 
most sensitive to the impurity concentration. For c — ► 1 
all curves are expected to join at a single value corre- 
sponding to the staggered magnetization of the planar 
Heisenberg AFM. For c = and J > 1.6 > J c , M q 
has been forced to zero in fig. QJb). This neglects finite 
size effects, namely that Mq 7^ 0, albeit small, even for 
c = since N = 24 is finite. In principle this value of 
Mq (J > J c ,c~ 0) can be read off by extrapolating to 
c — > for each J > J c in fig. Efb) . 

Finally we discuss the Chakravarty-Halperin-Nelson 
ratio 
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Figure 6: (a): Chakravarty-Halperin-Nelson ratio R of the 
total staggered structure factor and susceptibility vs. tem- 
perature at J = 2 > J c for c = 0, 0.03, 0.1, and 0.2, corre- 
sponding to dashed, solid, dashed-dotted and dashed-double- 
dotted. Result for c = in panel (a) from ref A Tempera- 
tures are 1/1024 < T < 1 (0.05 < T < 10) for c / (c = 0). 
Dotted line: unity, (b) and (c): denominator x (Q) an d nu- 
merator S (Q) from eqn. for R vs. temperature with c as 
in fig. Hj](a). Size of statistical errors is given by vertical bars 
in all panels. 



This ratio relates our analysis to that of the AFM non- 
linear ex-model (NLcrM) 1 ! 2 !^ . Regarding clean systems 
with a QCP between AFM-LRO and dimerization, this 
ratio has been studied in detail for the bilayer Heisen- 
berg and for the SKN21. In agreement with the 

NLcrM it was found that R = 1 in the (classical) high- 
T regime, as well as in AMF-LRO (renormalized classi- 
cal) regime for T — > 0. At the QCP the NLcrM requires 
R(T — > 0) ~ 1.09, which is also consistent with SSE at 
J = J c . While - to our knowledge - rigorous results are 
absent for a disorder induced AFM-LRO phase, as for 
J > J c , it is yet tempting to speculate that renormal- 
ized classical behavior will re-emerge. Since in the gaped 
phase of the clean system R diverges as T _1 this will lead 
to drastic variations in R versus c for J > J c . Indeed this 
is observed in fig. EJa), where R is shown versus T as c 
varies from to 0.2. Obviously R decreases markedly 
with increasing c and almost reaches 1 for c > 0.1 and 
T ~ 0.01. As for c = the zero temperature limit of R 
will be very sensitive to finite size gaps^LSS. This may 
cause the low-T increase seen in this figure, i.e. we specu- 
late that for larger systems R — 1 as T — > for all c ^ 0. 
This has to clarified by future analysis. For completeness, 
fig. Elb) and (c) relate the cancellation of the factor T 
in the denominator of R at c ^ to the temperature 
dependence of \ (Q)- 



which, in addition to S (Q) includes information on the 
longitudinal staggered susceptibility 



x(Q) 



l/T 



dr (m z Q (t) m z Q ) 



(6) 
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